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The increasing demand of electric power energy and the presence of 
disturbances can be identified as the factors of voltage instability condition in 
a power system. A secure and reliable power system should be considered to 
ensure smooth delivery of electricity to the consumers. A power system may 
experience undesired event such as voltage instability condition leading to 
voltage collapse or cascading collapse if the system experiences lack of 
reactive power support. Thus, to avoid blackout and cascaded tripping, 
load shedding is the last resort to prevent a total damage. Under Voltage 
Load Shedding (UVLS) scheme is one of the possible methods which can be 
conducted by the power system operators to avoid the occurrence of voltage 
instability condition. This paper presents the intelligent based technique for 
under voltage load shedding in power transmission systems. In this study, 
a computational based technique is developed in solving problem related to 
UVLS. The integration between a known computational intelligence-based 
technique termed as Evolutionary Programming (EP) with the under-voltage 
load shedding algorithm has been able to maintain the system operated 
within the acceptable voltage limit. Loss and minimum voltage control as the 
objective function implemented on the IEEE 30-Bus Reliability Test System 
(RTS) managed to optimally identify the optimal location and sizing for the 
load shedding scheme. Results from the studies, clearly indicate the 
feasibility of EP for load shedding scheme in loss and minimum voltage 
control in power system. 
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1. INTRODUCTION 


Today’s electric power systems are often subject to stress due to heavy loading conditions. 
In addition, the voltage instability is generally created due to component outage or sudden load increments. 
The current power system network has become more complex and heavily loaded due to growth in electricity 
demand, achieving maximum economic benefits and maximum transmission capacity efficiency. 
From literature it is observed that major power blackouts that have occurred around the world were caused by 
voltage instability. The increasing stress on the traditional power systems has resulted in increased 
blackouts [1]. One of the main causes of power blackouts is voltage instability which is attributed to 
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insufficient generation as well as transmission capacities. In situations where the transmission capacity is not 
enough for the system to operate normally, power systems employ load shedding mechanisms to reduce the 
load on the transmission network. Other major challenges faced by power system operators include a change 
in the nature of loads, the performance of the on-load tap changer transformer, the dependency on generation 
positioned remotely away from load centers, natural load growth, and the influence of protection and control 
systems. In these situations, the system may start moving toward collapse. However, before the voltage 
collapse point, the system enters into a contingency state. A system is said to be in a contingency state if the 
load on the system exceeds to a certain limit or a sudden change in the load value occurs leading to an outage 
of a transmission line or a generator [2]. 

The stability of a system is dependent on the voltage profile of the power network [3, 4]. When a 
system starts moving toward collapse, one of the indicators is the decline in bus voltages, as the voltage 
magnitude depends on the reactive power circulating in the system. Therefore, usually when the voltage 
magnitudes on various buses start to decline, motors are used to recover the voltage magnitude by increasing 
the amount of reactive power. However, in case of contingency conditions, the existing reactive power 
sources are not sufficient to stabilize the voltage profile. Moreover, factors such as unexpected load 
increments or component outages can worsen the situation further. Similarly, voltage instability may not only 
influence the local load area but may also spread to the adjacent areas in an interconnected power system, 
commonly known as cascading failures. Traditionally under-voltage load shedding (UVLS) [3] has been used 
to recover systems in contingency states and avoid a voltage collapse or blackout. UVLS has been proven to 
be a robust tool in stabilizing systems suffering from low voltage magnitudes [5, 6]. The success of UVLS in 
stabilizing a system depends on the optimality of the amount, time and location for load shedding. Shedding 
lesser or more than the required amount of load does not arrest voltage instability and may even lead to a 
voltage collapse or over frequency problems, respectively. Similarly, shedding load at the wrong place may 
cause unnecessary interruption, loss of customer trust, and the utility revenue [7]. The time instant at which 
load shedding needs to be performed is also very crucial as discussed in [8]. Some of the reliable 
optimization techniques for UVLS in smart grids are genetic algorithm (GA) and particle swarm optimization 
(PSO) as reported in [9, 23, 24]. To reduce the computational time, fast search technique can be adopted 
which involves optimization process. Several optimization techniques can be adopted in order to solve this 
problem such as; simulated annealing, linear programming, Genetic Algorithm (GA) and Evolutionary 
Programming (EP). GA and EP are two optimization techniques based on natural generation and specified 
under subdivisions of Evolutionary Computation (EC) in the Artificial Intelligence (AI) hierarchy. Two other 
methods under the EC that can be possibly adopted are Evolutionary Strategy (ES) and Genetic Programming 
(GP). The main difference in these techniques is normally in terms of the evolution methods. 

The proposed technique uses the bus voltage as objective function to select and shed the optimal 
load on the weak buses. The voltage stability margin is highly influenced by the weakest buses in an 
interconnected power system [11, 14, 16, 18]. Therefore, the identification of weak buses is necessary for the 
planning and operation of power systems. A weak bus is defined as a bus whose load-bearing capability is 
very low and the voltage magnitude is close to collapse identified by a low voltage value of the bus 
connected to it. 

This paper focuses on developing an intelligent algorithm for undervoltage load shedding, integrated 
with Evolutionary Programming. The proposed technique has been implemented on the IEEE 30-Bus RTS, 
resulting convincing results. 


2. RESEARCH METHOD 
2.1. Overview of Load Shedding 

Load shedding was identified as a corrective action against voltage collapse as proposed in [8]. 
In this work, voltage modal analysis provided a measure of the relative participation of each bus in voltage 
collapse scenario. The proposed EP technique revealed that load shedding could be initiated at the most 
critical locations, while the magnitude of the load to be shed could be determined from the reactive power 
deficiency and the sensitivity of deficiency to MW load provided by the interior point optimal power flow. 


2.2. Problem Formulation 

In this section, UVLS will represent as an optimization problem. Then this optimization problem is 
used in the algorithm for UVLS scheme as described in Section 2.3. The main objective for UVLS is to 
maximize the minimum voltage in the whole system while the sizing of reactive power to be shed and 
location will be determined by EP. The following equations present the equality and inequality constraints 
equations involved in this study. 
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As shown in (1) and (2) represent the equality constraints for active and reactive powers. On the 


other hand, (3) is the inequality constraint to maintain the voltage between the allowable range of minimum 
and maximum voltages. 


2.3. Proposed Method 


As mentioned in the previous section, EP involved initialization, fitness calculation, statistics, 
mutation, combination and selection. Detailed explanations on each step are described below. The overall EP 
process is represented in the form of flow chart as illustrated in Figure 1. 
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Figure 1. Flowchart of overall EP structure 


Step 1: Initialization: Generate an initial random number of variable x), Xo, ..... . Xx. These random numbers 
are the control variables which represent the locations and sizing of the load to be shed. 


Step 2. Fitness 1: Equation to be optimized. In this study, minimum voltage was used as the objective 
function and determined using the load flow program. 


Step 3: Mutation: Mutation was performed on the generated random numbers, x; to produce the 
offspring/children. 
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Step 4: Fitness 2: Calculation of fitness was conducted utilizing the mutated individuals or offsprings. 

Step 5: Combination: To combine parents and offspring in in cascade form which in turns doubled the 
number of individuals in the combined population. 

Step 6: Selection: To identify the candidates to be transcribed into the next generation. 

Step 7: Convergence Test: To determine the stopping criteria of the evolution. 


3. RESULTS AND ANALYSIS 

In this section, an optimization engine and UVLS algorithm were developed in integrated form in 
MATLAB. The validation processes were conducted on the IEEE 30-Bus RTS as shown in the single-line 
diagram in Figure 2. In this study, 4 cases were considered as described in the succeeding sections. 








Figure 2. IEEE 30-Bus Reliability Test System (RTS) 


3.1. Case 1: Single Bus Load Shedding 

In this case, six buses have been considered to undergo load variation independently. Buses 14, 15, 
16, 18, 21 and 24 are the buses subjected to 30 MVAR load at a time. One load bus was optimally identified 
by EP to experience the load shedding. 


Table 1. Load Shedding Under Single Load 








Tested Bus Voltage before load Voltage after load Bus Amount of load 
(Qi=30MVAR) shedding (p.u.) shedding (p.u.) Location shed (MVAR) 

14 0.9740 0.9901 20 16.6913 

15 0.9807 1.0035 25 19.8207 

16 0.9825 1.0016 25 19.8207 

18 0.9487 0.9837 19 15.3299 

21 0.9825 1.0037 25 19.8207 

24 0.9643 1.0006 25 19.8207 





Table 1 tabulates the results of shedding | load with the location and amount of reactive power to be 
shed. Among the 6 buses, bus 18 exhibits the lowest bus voltage 0.9487 p.u. before the implementation of 
EP-UVLS scheme. This voltage is improved to-0.9837 p.u. with the implementation of UVLS. The proposed 
EP has identified bus 19 to undergo 15.3299 MVAR reactive load to be shed to improve this voltage value. 
For Qa4=30 MVAR, the implementation of integrated EP-UVLS managed to improve the voltage value from 
0.9643 p.u. to 1.0006 p.u., which requires an optimal value of 19.8207 MVAR to be shed at bus 25 as the 
corresponding load bus. Looking at the whole table, it shows that bus 25 is identified as the most frequent 
load bus to experience load shed regardless of tested bus which was subjected to 30 MVAR load. 
This information could be useful to power system operators responsible for their UVLS scheme. On the other 
hand, Figure 3 shows another presentation of results to demonstrate the difference in profile before and after 
EP-UVLS implementation. Apparently, all the participated buses experience voltage improvement with EP- 
ULVS scheme. Figure 3 presents the difference of voltage before the pre and post-load shedding for all the 
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participating buses in the system. In general, all the voltage values at all the participating load buses have 
been increased with the implementation of load shedding scheme. 
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Figure 3. Voltage improvement for single load shedding 


3.2. Case 2: 2 Buses Load Shedding 
In this case, the same 6 load buses i.e. buses 14, 15, 16, 18, 21 and 24 were independently subjected 


to 30 MVAR. Table 2 tabulates the results for load shedding at 2 buses. In general, all the voltage values 
have been improved with 2 load buses to be shed. The voltage at bus 18 has been improved within acceptable 


operating limit even though the voltage is still less than 1.0 p.u.. 


The proposed EP-UVLS algorithm had indicated that bus 17 and bus 25 as the frequent locations to 
experience load shedding for all the tested buses. The amount of reactive power to be shed are 12.6528 
MVAR and 19.8207 MVAR for buses 17 and 25 respectively as shown in the table. The results for voltage 
increment at other tested buses can be referred to the same table. Figure 4 presents the difference of voltage 
before the pre and post-load shedding for all the participating buses in the system. In general, all the voltage 
values at all the participating load buses have been increased with the implementation of load 


shedding scheme. 


Table 2. Load Shedding Involving 2 Load Buses 











Voltage before Amount of Amount of 
Roeteeen load shedding V0Wage alter load Bus Toadshed BUS toad shed 
= (pu) g (p.u.) ocation (MVAR) ocation (MVAR) 
14 0.974 1.0034 24 17.3649 23 14.7245 
15 0.9807 1.0043 17 12.6528 25 19.8207 
16 0.9825 1.0042 17 12.6528 25 19.8207 
18 0.9487 0.9883 19 18.8196 11 15.2303 
21 0.9825 1.0044 17 12.6528 25 19.8207 
24 0.9643 1.0041 17 12.6528 25 19.8207 
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Figure 4. Voltage improvement for load shedding involving 2 load buses 
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3.3. Case 3: 3 Buses Load Shedding 

In this case, the same buses (buses 14, 15, 16, 18, 21 and 24) were tested for load shedding scheme 
involving 3 load buses. In Table 3, bus 18 which experiences the lowest voltage value has witnessed voltage 
increment from 0.9487 p.u. to 0.9964 p.u. with load shedding involving 3 load buses. For this load bus, buses 
27, 18 and 19 are the optimal location identified by the EP-UVLS scheme; with the corresponding load to be 
shed worth 5.3863 MVAR, 10.7461 MVAR and 5.74420 MVAR respectively. This is a good indication since 
the system can still operate within the secure voltage limit even only small reactive loads to be shed at the 
corresponding buses. Apparently, buses 20, 25 and 28 are the frequent buses recorded as the buses to be shed 
in this system. One important observation in this case is that, buses 15, 16, 21 and 24 witness identical 
optimal locations and sizings for load shedding as indicated in Table 3. 

Bus 14 almost reaches its maximum value as the table shows the proposed EP gives the same 
location for the 2" and the 3™ locations. The variation of voltage value for all the participating load buses for 
three load shedding scheme is presented in Figure 5. 


Table 3. Load Shedding Involving 3 Load Buses 








Voltage Amount Amount Amount 
Voltage before after load Bus of load Bus of load Bus of load 
Tested Bus load shedding shedding _ Location shed Location shed Location shed 
(Qd=30MVAR) (p.u.) (p.u.) 1 (MVAR) 2 (MVAR) 3 (MVAR) 
14 0.974 1.0018 23 14.7245 25 19.9892 28 4.6639 
15 0.9807 1.0051 20 19.1072 25 11.0917 28 12.4970 
16 0.9825 1.005 20 19.1072 25 11.0917 28 12.4970 
18 0.9487 0.9964 27 5.3863 18 10.7461 19 5.74420 
21 0.9825 1.0052 20 19.1072 25 11.0917 28 12.4970 
24 0.9643 1.0049 20 19.1072 25 11.0917 28 12.4970 
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Figure 5. Voltage improvement for load shedding involving 3 load buses 


3.4. Case 4: 4 Buses Load Shedding 

For the last case, 4 buses load shedding scheme involving the similar load buses as the tested buses 
as those in Case 1, Case 2 and Case 3. The reason of conducting several cases involving similar tested load 
buses is to see the consistency in terms of voltage increment and sizing of loads. The results for voltage 
increment, optimal location and sizing are tabulated in Table 4. The lowest voltage experienced by the 
system when bus 18 was subjected to 30 MVAR has been significantly improved to 0.9994 p.u. when 4 buses 
are shed. These buses are buses 14, 21, 11 and 20. The corresponding reactive loads to be shed are 3.0744 
MVAR, 16.0481 MVAR, 10.689 MVAR and 19.1072 MVAR. Surprisingly bus 21 requires significantly the 
largest reactive load to be shed worth 16.0481 MVAR as compared other loads. However, the first bus, i.e. 
bus 14 only requires small reactive load to be shed. In general, it is acceptable since the proposed EP-ULVS 
is intelligent enough to identify the location and sizing of the reactive load to be shed. Apparently, buses 25, 
13, 24 and 23 are the most frequent buses to experience load shedding. This will become a good indicator for 
the system when all the 6 buses are involved. 
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Table 4. Load Shedding Involving 4 Load Buses 








Tested Bus Voltage Voltage Amount Amount Amount Amount 
(Qy=30MVAR before load after load Bus of load Bus of load Bus of load Bus of load 
shedding shedding Loc 1 shed Loc 2 shed Loc 3 shed Loc 4 shed 

) (p.u.) (p.u.) (MVAR) (MVAR) (MVAR) (MVAR) 

14 0.974 1.0051 25 16.5563 13 0.3174 24 17.3649 23 14.7245 

15 0.9807 1.0051 25 16.5563 13 0.3174 24 17.3649 23 14.7245 

16 0.9825 1.005 25 16.5563 13 0.3174 24 17.3649 23 14.7245 

18 0.9487 0.9994 14 3.0744 21 16.0481 11 10.689 20 19.1072 

21 0.9643 1.0052 25 16.5563 13 0.3174 24 17.3649 23 14.7245 

24 0.9643 1.0050 25 16.5563 13 0.3174 24 17.3649 23 14.7245 





Further results presentation for this case is depicted in Figure 6. This figure demonsrates that the 
implementation of load shedding involving 4 buses has managed to improve the voltage profile at all the 
participating buses in the system. The idea is to show the graphical presentation for voltage increment at all 
the participating load buses to be tested. In general, all the voltage levels are within the operable voltage limit 
Le. 0.95<V,,<1.05 p.u.as stipulated by most standards in power utilities. 
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Figure 6. Voltage improvement for load shedding involving 4 load buses 


4. CONCLUSION 

This paper has presented intelligent based technique for under voltage load shedding in power 
transmission systems. In this study, an integrated intelligent which combined EP and UVLS termed as EP- 
UVLS has been developed and validated on the IEEE 30-Bus Reliability Test System (RTS). The proposed 
EP-UVLS algorithm has significantly improved the voltage profile of a power system which is feasible to 
avoid a power blackout. The proposed integrated EP-UVLS technique has shown its capability to find the 
optimal location and sizes of reactive power (MVAR) with the objective to improve the voltage profile in the 
transmission system. All of the objectives in this study have been significantly achieved. However, there is 
still room for improvement. For instance, new technique such as Monte-Carlo can be integrated in the 
initialization process in order to generate the random number. The proposed EP-UVLS can also be improved 
in terms of its computational time by integrating fast search technique subject to further study. Furthermore, 
integration of renewable energy (RE) can also be installed in the system to see the impact of its presence in 
load shedding scheme. Results from the study can also be taken as early indication and information for the 
planning in power system. 
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